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Propagation of the photoelectron wave packet in an attosecond streaking experiment 
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Laser-assisted photoemission from a solid is considered within a numerically exactly solvable one¬ 
dimensional model of a crystal. The effect of the inelastic scattering and of the finite duration of 
the pump pulse on the photoelectron dynamics is elucidated. The phenomenological result that 
the photoexcited wave packet moves with the group velocity dE/dk and traverses on average the 
distance equal to the mean free path is found to hold for energies far from the spectral gaps of the 
hnal state band structure. On the contrary, close to a spectral gap the photoelectron is found to 
move during the excitation by the pump pulse with a velocity higher than the group velocity. 


INTRODUCTION 


Transport properties of electron wave packets underlie 
the functioning of electronic devices and are an impor¬ 
tant factor in photoemission spectroscopies and electron 
diffraction techniques. In a solid, as opposed to a sin¬ 
gle atom, special role is played by inelastic scattering 
giving rise to a limited electron lifetime and mean free 
path. The inelastic scattering is phenomenologically de¬ 
scribed as electron absorption [ij , and its iimlications for 
stationary processes are well understood A much 

less studied question is how the wave packet evolves and 
propagates on a time scale comparable to its lifetime. 
Experimental access to such ultrafast processes is given 
by modern time-resolved photoelectron spectroscopies, in 
particular by the streaking method |8l4l2l|. in which a 
subfemtosecond time resolution is achieved by mapping 
time to energy using a strong laser field: the electron 
wave packet created by an ultrashort pulse of extreme 
ultraviolet radiation (XUV) is accelerated by the super¬ 
imposed laser field, and the energy by which its spectrum 
shifts up or down depends on the electron release time tx 
relative to the temporal profile of the laser field 

Various models have been developed that treated the 
problem W the time-dependent Schrddinger equation 
(TDSE) or by applying a classical model [l^ . 

The majority of the studies consider a solid in which the 
low-frequency laser field is strongly damped by the di¬ 
electric response [l^ , so the electron needs to travel some 
distance before it gets exposed to the streaking field. 


MEAN FREE PATH OR LIFETIME? 


Assuming an instantaneous XUV excitation, the time 
of flight from the depth where the photoelectron is cre¬ 


ated to the surface is determined by the group velocity of 
the wave packet. Because of the inelastic scattering the 
photoelectrons excited deeper in the crystal have smaller 
probability to reach the detector, which is assumed to 
scale with the depth z as exp(z/A). Then the average 
depth traversed by the electrons is just the mean free 
path (MFP) A, and the time to reach the surface can be 
estimated as r = X/v, where v is the group velocity of 
the emitted wave packet [l2| . Often it is sufficient to use 
the values of A averaged over an energy interval of several 
eVs that can be inferred from the well-known universal 



FIG. 1: Dependence of optical potential U on energy. Small 
circles are the values empirically derived from ab initio analy¬ 
sis of ARPES spectra (full circles) and VLEED spectra (open 
circles). Large circles are the ab initio GW calculation for 
W(IOO) (FLEUR-SPEX). The dependence V/E) is well fitted 
by a function 7 arctan(i5//3) with /? = 80 eV and 7 = 5.3 eV 
for W(IOO), 7 = 4 eV for Al(lll), and 7 = 6.2 eV for TiTe 2 . 
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X{E) curve. This is, however, not sufficient for ARPES 
because MFP may strongly vary over a few eV interval 
as a result of the specific band structure. At the same 
time, experience with ARPES and VLEED (very low en¬ 
ergy electron diffraction) suggests that the photo¬ 

electron final-state lifetime changes rather smoothly with 
energy, see Fig.[T](a). Thus, it is reasonable to rely on the 
average values for the lifetime, and to derive the rapidly 
varying MFP as X = vt = u/V), where Vi is the inverse 
lifetime (optical potential), which enters the Hamilto¬ 
nian through the imaginary term —iVi. For a nearly free 
electron, the wave vector k acquires an imaginary part 
K = mVi/kh^, which yields A = 1/2 k, and the time to 
reach the surface is then 



Note that the group velocity cancels, so r depends on 
energy through the function Vi{E). Vast empirical data 
suggest that the Vi{E) curves are similar for different 
materials: Vi grows with energy with similar rate, see 
Fig. Ha). Because the empirically derived data may be 
affected by various extrinsic factors and by the imper¬ 
fection of the apparatus it is important to theoretically 
corroborate this observation. Microscopically, optical po¬ 
tential is associated with the imaginary part of the self¬ 
energy. Figure [TJa) shows an ab initio GW calculation 
for a high-energy conducting band of W(llO). It is seen 
to agree well with the empirical results, which supports 
our understanding of the electron absorption as coming 
predominantly from the electron-electron interaction and 
corroborates the empirical estimate of its growth with en¬ 
ergy. 

Thus, within this simplest approach the arrival delay of 
photoelectrons with final energy E relative to those with 
energy E -\- e behaves roughly as At{E) ~ e/E{E — e), 
see Fig. HF): the electron with a lower energy comes 
later, and the delay At decreases as the final state en¬ 
ergy grows. The recent measurements by Neppl [13 on 
W(llO) of the phase shift between the spectrograms of 
valence 5d and semi-core 4/ states qualitatively follow 
this trend [open circles in Fig. [TJb)], but the absolute 
values in the experiment are about three times larger. 
At the same time, the value of 5 ± 20 as measured on 
Mg(OOOl) at hui = 118 eV is 4 times smaller than ex¬ 
pected from the simple MFP-based theory. 

FUNDAMENTAL QUESTIONS 

Such drastic disagreement in absolute values calls for 
the search of mechanisms that, in addition to the finite 
lifetime, affect the transport of the photoelectron to the 
surface. A number of fundamental questions arise: What 
is the role of the finite duration of the pump pulse, i.e., at 
what point does the wave packet start its motion with the 


group velocity? Does the electron have the group veloc¬ 
ity dictated by the band structure if it travels only the 
distance of a few atomic layers? How does the packet 
propagate when its energy falls into a forbidden gap, 
where the group velocity cannot be defined as dE/dkl 
Can the electron excited at the outermost atomic plane 
be thought of as not feeling the band structure because 
it moves away from the surface, towards vacuum? 

Apart from that, an important theoretical question 
arises of whether the absorbing potential —iV can be 
used in the time-dependent Schrodinger equation to ade¬ 
quately describe both the photocurrent attenuation and 
the implications for the transport: inelastic scattering is 
a stochastic process, so the photocurrent is an incoher¬ 
ent overlap of randomly dephased contributions, whereas 
the optical-potential-damped wave packet remains fully 
coherent, which may have implications both for its trans¬ 
port and for its acceleration by the laser field. 

THE MODEL 

To resolve these questions, we need a model that does 
not incorporate any assumptions about the answers and, 
at the same time, can be treated numerically exactly. 
For the present proof-of-concept calculation we use a 
one-dimensional model of a crystal, which treats pho¬ 
toexcitation, transport, and streaking fully quantum- 
mechanically without resorting to a perturbational treat¬ 
ment of any of the terms in the Hamiltonian or to a 
heuristic separation of the XUV and the laser pulses. We 
employ a particle-in-the-box method: the box comprises 
a thick crystal slab lying on a structureless substrate [the 
piecewise constant potential V{z) in Fig. Hs-)] fFe 
vacuum half-space. In Ref. [2^ this model was applied to 
photoelectron streaking by a spatially uniform laser field 
in the absence of inelastic scattering. The present model 
allows for a rapid decay of the laser field into the solid 
and includes the inelastic scattering in two alternative 
ways: via the optical potential and by a straightforward 
ensemble averaging over random configurations. 

THE PROCEDURE 

We perform a series of numerical experiments, in which 
the system is excited by an XUV pulse of duration 
Dx = 1 fs and is simultaneously acted upon by the 
laser pulse of duration Dl = 5 fs. Fig. HF)- Both 
light fields are given by the dipole operator z, and the 
laser field is screened by multiplying it by a smooth 
step function 9{z). The temporal envelopes of both 
pulses are of the form cos^ (nt/D). The laser pulse is 
Shit) = sinD(t—A) cos^[7r(t—A)/!?!,], with photon en¬ 
ergy D = 1.65 eV and amplitude £■“ = 2 x 10^ V/cm. The 
TDSE is solved with the split-operator technique in ma- 
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FIG. 2: (a) Crystal potential V{z) with a defect at the third layer, (b) Superposition of the XUV and the laser pulses, (c) 

Streaked spectra for Tiui = 78 eV with absorption (lower curves) and with random collisions (upper curves) for II — tx = 1600 as 
(solid lines) and 200 as (dashed lines), (d) Band structure with the periodic potential V{z) and energy dependence of the 
cross-section of emission from a localized state at a defect, (e) and (f) Streaking curves for huj = 78 and 90 eV with random 
collisions (e) and with an absorption (f). (g)-(i) Surface arrival time as function of the final state energy. Small circles connected 
by solid lines in graph (g) and dashed lines in (h) and (i) are streaking results with absorption. Large circles in graphs (g) and 
(i) are streaking with random collisions. Solid lines in graphs (h) and (i) are from the equation of motion in vacuum. 


trix form in terms of exact eigenfunctions (discrete spec¬ 
trum) of the unperturbed Hamiltonian H = fp' j2m -\- V, 
so the crystal potential is fully taken into account for 
both initial and final states. The spectrum of H is trun¬ 
cated at 200 eV above the vacuum level, which ensures 
the convergence of all the results 221. 


Apart from Hq and the two external fields the Hamilto¬ 
nian includes inelastic scattering, which is either a static 
absorbing potential —iVi[l — 9{z)\ or a stochastic per¬ 
turbation ( 7 (t)s„(z)[l — 0(z)], where q{t) and Sniz) are 
random functions of t and z, respectively. Here q(t) 
is fixed and is the nth sample in a random ensem¬ 
ble. The ultimate spectrum is then the average over A^r 
random walks, Fig. He). The peak displacement AE 
from its laser-free position as a function of the time shift 
At = E — tx between the XUV and the laser pulse gives 
the key to the temporal information [^ : by fitting the 
measured AE(At) points with the momentum transfer 

OO 

function p{t) = ejdt£i^(t — tx), we infer the arrival time 

T 

of the electron in vacuum r, see Fig. H®) the motion 
in the random potential and Fig. [5](f) for the motion in 
the absorbing potential. 


RESULT: XUV PULSE SHIETS PACKET 


For a detailed study of the photoelectron transport to 
the surface it is convenient to exactly know its initial po¬ 
sition in space. This is achieved by introducing a small 
defect at one of the layers and considering the photoe¬ 
mission from the localized state at the defect. 

Figures me) and mf) show streaking curves for the 
initial state localized at the third layer for huj = 78 and 
90 eV. Note that the results by the phenomenological 
absorbing potential (f) and by the microscopic random 
collisions (e) perfectly agree. Counterintuitively, the elec¬ 
tron at the higher energy arrives 40 as later than at the 
lower energy: r = 150 as at 90 eV and r = 111 as at 
78 eV. Figure ms) shows the arrival time of the electron 
initially at the 2nd, 3rd, and 4th layer for photon energies 
from 65 to 159 eV. The most striking are the two minima 
at huj = 80 and at 143 eV {E = 39 and at 102 eV), at 
which the arrival time even shows negative values. 

The minima B and C are located at the lower edges 
of the band gaps, and minimum A occurs at the energy 
where the ionization cross-section vanishes, see Fig.^d). 
There the t{E) curve strongly deviates from that derived 
from the group velocity, t{E) = d/v{E) [dashed curve in 
Fig. mg)l- To prove that the discrepancy is not due to 
the rather indirect streaking method to measure r, we 
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FIG. 3: Temporal evolution of the photoelectron spectrum 
from a localized state at hto — 80 eV (left) and Tiui = 100 eV 
(right). The black bars at 46 and 50 eV indicate a spectral 
gap, and the red bar at 39 eV indicates the Cooper minimum. 


approach the dephasing of the wave packet is neglected, 
and the question arises if this may bring about any ar¬ 
tifacts. This simplification turns out to be not crucial: 
Figure [2Dg) shows the t{E) data by statistical averaging 
over 144 configurations (large circles) for emission from 
the 2nd and 4th layer and Fig. Hi) for the 1st layer. The 
parameters of the random perturbation are chosen such 
that it gives the same MFP as V) = 1.7 eV. The dephas¬ 
ing is seen to cause a broadening of the t{E) features, 
but otherwise the behavior of the incoherent ensemble is 
in perfect accord with that of the coherent packet mov¬ 
ing in the absorbing medium. That such two very dif¬ 
ferent systems agree to within fine details suggests that 
this advancement phenomenon is robust, and it should 
be expected to occur in real solids. 


CONCLUSIONS 


employ an alternative method: we switch off the laser 
and follow the propagation of the wave packet up to a 
distance of 500 a.u. away from the surface. By measur¬ 
ing the center of gravity of the packet in a number of 
time points we determine its equation of motion in vac¬ 
uum Zq + vt and obtain the time point to at which it has 
arrived at the surface. The to{E) and t{E) curves are 
seen to agree well. Fig. [2])h), which means that the clock 
provided by the streaking technique is accurate, and that 
the fast delivery of the photoelectron to the surface is ef¬ 
fected during the excitation process. Figure [3] shows the 
time evolution of the spectrum at the t(E) minimum, 
huj = 80 eV, and at huj = 100 eV, where t{E) well agrees 
with the instantaneous approximation d/v{E). The spec¬ 
tral evolution maps in the two cases are qualitatively dif¬ 
ferent: while at 100 eV the spectrum rapidly concentrates 
around E = huj + Ei, at 80 eV the evolution is much 
slower: at t = 200 as the intensity is still spread over a 
40 eV wide spectral interval. Clearly, while the spectral 
coefficients il){E,t) of the packet JdE ilj{E,t) \ E ) keep 
changing, the velocity of its gravity center may deviate 
from the weighted group velocity of the stationary waves 
\ E ), and during the XUV pulse the wave packet may 
travel a distance comparable with the mean free path. 

Such unusual behavior happens every time when the 
central energy of the wave packet approaches an intensity 
minimum, be it a vanishing matrix element (Cooper min¬ 
imum) or a spectral gap. For the emission from the first 
layer, at its Cooper minimum {Huj = 77 eV) the XUV 
pulse, on the contrary, detains the packet in the solid 
causing a delay of ^ 200 as. All the other features of the 
t(E) curve are similar to those of the emission from the 
deeper layers. 

Figure shows that the inelastic scattering plays 
an important role in the formation of the spatial shape 
and transport of the packet: the displacement effect is 
more pronounced at larger V)- In this phenomenological 


The present calculation has revealed a peculiar dynam¬ 
ics of photoelectrons excited by an ultra-short pulse. The 
interesting features arise from the interplay between the 
elastic scattering from the crystal lattice and inelastic 
scattering that causes the damping of the electron wave. 
Three aspects can be distinguished: first, the effect of 
inelastic scattering is to increase the average velocity of 
the wave packet, which is natural because it is the slower 
components that are damped stronger, as they spend 
more time in the solid. The range of the advancement 
effect depends on the interaction with the crystal lattice: 
at the special points, and around the gaps it may exceed 
200 as, whereas in the free-electron-like region it reduces 
to 20 as, which is still comparable with the measured rel¬ 
ative delays, see Fig.[TDb). This effect can be understood 
in terms of the group velocity of Bloch electrons. 

At certain photon energies, however, owing to a 
complicated non-free-electron-like structure of the final 
states, the wave packet performs a complicated motion, 
and it may approach the surface much faster than an 
instantaneously created packet. This may also lead to 
an interesting effect that the electron starting from the 
outermost layer is overtaken by the electron coming from 
the depth of the crystal. 

Close to the spectral gaps the temporal shift of the 
spectrogram does not necessarily agree with the equation 
of motion of the center of gravity of the wave packet, 
but apart from the special points discussed above the 
discrepancy of the two methods does not exceed 20 as. 
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